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ABSTRACT 


0-C diagram is a useful technique to analyse the period changes of a pulsator 
by using the maximum (or minimum) value points that have been obtained from 
the historical data. But if an object is a double-mode or multi-mode pulsator, 
the extreme value points are the results of all the modes other than just the 
fundamental mode. We will discuss these situations and give out some criteria 
to judge whether 0-C diagram can be used under these situations. 

Subject headings: techniques: 0-C diagram —stars: variables: double-mode pulsators 
— stars: variables: multi-mode pulsators 
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1. Introduction 


0-C means 0[bserved] minus C[alculated], it is a useful method that we always use to 
analyse the period change of variables. More over, we do not plan to discuss the details how 


to do 0-C diagram, which would be found in 


SterkenI (1200 


3 ). 


The usage of 0-C diagram in variables depends on the maximum value points of 
the light curves we observed, but these points are always unstable when the object is a 
double-mode or multi-mode pulsator. All the modes coupling together and the maximum 
points are the comprehensive effects of them all. Of course, we can use the other methods 
to analyse the period change of the pulsator for example the Fourier analysis for each data 
set. But in most cases, these methods need a lot of high quality data, which is always 
difficult to obtain. Furthermore, in the historical observational data of most variables, the 
maximum value points are the most credible and always the only credible parts. We should 
make full use of the data. As a result, we give out some criteria to decide whether the 
classical 0-C diagram can be used to these pulsators. 


In the following sections, we present a detailed discussion of how to estimate 
transformation of the maximum value points and give some criteria which can be compared 
with the observational errors and fit errors of 0-C diagram. We also assume that a global 
pulsation analysis has already been made, and the frequencies of the pulsator have been 
obtained. 


The organization of the paper is: Section 2 analyses the transformation of the maximum 
values and gives out the criteria; Section 3 applies the criteria to a specific example; Section 
4 gives our conclusions. 


^Send offprint request to: jnfu@bnu.edu.cn 
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2. The error estimation 

2.1. Combination of two different sinusoidal functions 


Let, 

1. fo{t) = Ao cos{ujot + 00); 

2. fi{t) = Aicos(a;it+ 0i); 

3- /o,l(^) ~ ^0 COS^Uot + 0o) 4“ COS^OJit + 0i). 

In the above equations, we choose /o and /i satisfying < 1. In most conditions, 
Aq > Ai and Uq < Ui, but this is not always the case. 

The main problem can be reduced to the problem below: let 

. = 0 and < 0; 

• /o,i(^o)i) = 0 and < 0, 

and then we want to get and which are the maximum value points of /o and 
/o,i, respectively. 

For any value of t™, we choose the that is closest to t™, and let 5T'^ = — 1'^]\ 

And we are here to want to estimate AT*” = max{5T^}. 


2.1.1. coi = nuJo,n ^ 

In this case, fo,i{t) is the combination of the harmonics of cuq, which do not influence 
the distance between two maximum values. Thus: (JT™ = 0. 
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2 . 1 . 2 . cui = aa;o,a e M+\ N+ 

The primary problem can be reduced to the problem below: 

1. fo{t) = Aocos{ujot)] 

2. fi{t) = cos(a;it + ^i); 

3. = ^0 cos(a;ot) + cos(a;it + 0i); 

the left are the same as before. All the phase effects attribute to 0i, 0i G M. Note that 
when 01 is a free parameter, all the other cycles can be represented by the cycle near t = 0 
if we choose a proper 0 i. 

We now focus on the maximum value near t = 0. We know that = 0, and let us 
estimate 

Let /o i(t) = 0, leading to: 

AqUq sin(ci;ot) + AiUi sin(ci;it + 0i) = 0. (1) 

Then, 

sin(a;ot) = —+ 0i) < (2) 

AqUq AqOJq 

When the equal sign is satished, Aisin{ujit + 0i) reaches its maximum value, and 
sin(a;ot) = Since sinusoidal function is monotonous in [—vr, tt], f^i reaches its largest 

value now. At the same time, = 0, therefore we can get: 


AT”^ = max{5m = — arcsin(4^). 

Wo Aowo' 


( 3 ) 
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2.2. The transformation of two adjacent maximum values 

Because all the values of are around the corresponding values of and if oji/ojq 
is a irrational number, then all the values in [0, AT”*] will be reached in some cycles (Of 
course, ojq/uji ate always in this case). Averagely speaking, we have 

E„(lg(»)-tgi(" + l)l) ^ S„(IC(n)-C(n + l)l) _ 27 r^ 

n^oo n n^oo n CJq 

In the equation, and t™(n) imply the n-th maximum values of /o,i and /o from 

t = 0. 

When we use 0-C method practically, the transformation of two adjacent maximum 
values = ||t[J]'i(n) — + 1)| — ^|) are more important than the absolute 

transformation (AT”*) relative to In other words, the upper limit of distance between 
two adjacent maximum values At”* = max{(5t”*} should be estimated. 

In this case, let: 

• fo(t) = Aocos{uot), /o(t) = -Aoa;osin(a;ot); 

• flit) = Ai cos(a;it + 0i), /((t) = -AiUi sin(a;it + 0i). 

In general conditions, 0o and 0i are not important if we only want to get a upper limit 
of (5t”*. 

Let us calculate the next period after t = 0: 


/o(^ + ^o) — —AqUo sin(a;o(^ + ^)); 

/((t + Ti) = —AiiOi sin(a;i(t + ^) + 0i). 




OJO 
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For ojq < cui, ^ let A(/) = 

Then, 


A 

OJI ■ 


• AqUo sin(a;ot) = —AiUi sin(ci;i(t + A0) + 0i). 

This can be reduced to the condition in section 2.1, and the solution of this equation near 
t = 0 stands for the between the hrst and the second maximum values of /o,i. 

Obviously, < AT™" and changes along with A0 periodically. 

Note that AT”^ = ^ arcsin(^^), and A0 depends on the cuq and cui, we push out a 
estimation of (For every two adjacent cycles, the phase differences would be A0.): 

5t^ = AT- sin(^A0) < AT^. (5) 

As a result, Af— could be given as: 

At™ = maxi*"*} = Ar"sm(,^A(^ - h-Ar’"sm(-T = AT"sin{^A4>-^) + AT"'. 

( 6 ) 


2.3. Combination of three different sinusoidal functions 

2.3.1. Harmonics 


For harmonic terms, the maximum value points of cuq would not be transformed. 




2.3.2. Three absolute modes 


• Estimation of AT*” 

A proper estimation can be given as follows, if we use the same method as before. 


AqcJq sin^cugf) — —AlOUl sin^cj^f -|- 0i) — A 2 CU 2 sin(cj2t 02) A -|- A 2 CU 2 I (7) 


If 


AiCJi + A2CU2 
AqCUo 


< 1 , 


( 8 ) 


then we can get 


AT" = 

UJq AqUJq 


( 9 ) 


• Estimation of At™' 

If we take the total effects as the additional effects of uji and a ;2 independently, then 
we can use their linear superposition as the estimation of Af”^. Consequently, the 
estimation should be 


At- = ATf sin(^A0i - ^) + ATf + AT^ sin(^A02 -\) + AT^. (10) 


AT™, AT™, A01 and A02 are dehned similarly. 


2.4. Generalization 
2.4.1. Harmonics 

In practical usage, all the harmonic terms of ujp {p G N) could be estimated totally as 
one trigonometric function for convenience. We use one sine function to estimate the total 
effects of the all the terms of ujp. 

fp{t) = Ap cos{u:pt + 0p) fp{t) = ^ Ap^q cos{qUpt + (pp,q).q G N, g > 1 

9 


( 11 ) 
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In this equation, Ap = maxt(zmathbbR{fp{t)} and ensures fp{t) = fp{t) when t = 

(t™ is the maximum value points of fp{t)). 

Because we use this function to estimate the maximum values, and in most cases, 
Ap^q/Ap^q_i > 5, fpi--^) would be a good approximation of Ap and (pp^i would be a good 
approximation of (pp. 

As a consequent, we choose 


fpi^') ApCosi^ujpt (pp)- ( 12 ) 

to reprensent the effects of fp{t) (sometimes we can also choose Ap = Ylq-^p,q 
approximation of Ap). 


2.4-2. Multi-mode combinations 


We just consider a p — mode (p > 1) pulsator as an example. Using the method above, 
we obtain the generalization as follows. 

If 

dloa;o'" ’ 

then we can give another estimation; 

AT™ = — arcsin( V 4^) 

Wo AquJo 


(14) 


(13) 


The estimation of transformation of adjacent maximum values should be: 

At™ = - |) + AT”). (15) 

P ^ 

In the equation above, AT™ = ^ arcsin(^^) and Acpp = ^ 
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3. An example (AE UMa) 

AE UMa is a SX Phe star which has two modes of pulsation. Because of its 
double-mode pulsation, it is doubtful that whether the classical 0-C method can be used 
to analyse the period change of its fundamental mode. The key is that the maximum 
values are the results of two mode together, the hrst overtone and the beat mode cause 
the transformation of the maximum values. Now, let’s estimate the upper limit of the 
transformation by using the criteria introduced above. 


3.1. Estimation of and At^‘ 


First, using the data from iNiu et al.l (120131) . Table 3, we collect each pulsation terms 
cuo and oji. 

Let’s deal with the combinations of cuq and coi (Amojo.ntJi cos((ma;o + nui)t + 
m,n ^ Z). 


All these terms can be transform to the form 


-^muJo,nu}i COs((?7l(no T TlUJl)t -|- (j)mLOo,nu}i') 

sm{uit + ipmujo ) cos(a;ot -h (pmu^o )+ 

sin(a;ot + V^mojo ^nuji ) cos(a;it -t- ipmwo ^nuji ). (16) 

Ill the epuation, '~PmuiQ,nLj\ 2^4^rruxiQ,nw\ ■ 2 )' 

And we can add the term Araojo,nm sm{uit + (pma;o,naji) cos(a;ot (Pmujo,nuji) to fo(t) and 
the term Amcoo ,nLJi Sm{Uot -f ^Pmcoo ,nuji ) COs(^UJit ^mujo ,noJi ) to 

In addition, we only choose the frequencies whose amplitude larger than lOmma^f, 
since the smaller ones are considered to influence the result slightly. 
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The reduced /o and fi can be obtained by Ap = and (j)p^i. According to the 

calculation we introduced: AT™ = 0.002728, At™ = 0.001504. This result implies that the 
upper limit of the transformation between adjacent maximum values is in the same level 
with the observational errors and the £t error from 0-C diagram, which have the value of 
0.002083 and 0.00244 respectively. 

As a result, the classical 0-C method can be used to determine the period changes in 
this situation, but just for the fundamental mode. 


4. Conclusion 

In this paper, we give out the criteria which can tell one whether the classical 0-C 
method could be used to analyse the maximum values observed from a multi-mode pulsator. 
According to comparing the upper limit estimation of the maximum value transformation 
(AT™, At™) and the observational errors or £t errors of 0-C diagram, we can make a 
decision whether the 0-C method can be used effectively. If At™ < a, we think it is 
meaningful to use the 0-C diagram to do analysis. 
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